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Abstract. The cTTS TWHya has been observed with high-resolution X-ray spectrometers. Previously found high 
densities inferred from He-like f/i triplets strongly suggested the detected X-ray emission to be dominated by an 
accretion shock. Because of their radiation field dependence He-like f/i ratios do not provide unambiguous density 
diagnostics. Here we present additional evidence for high densities from ratios of Fexvn lines. Key Fexvn line 
ratios in TWHya deviate from theoretical expectations at low densities as well as from the same measurements 
in a large sample of stellar coronae. However, a quantitative assessment of densities is difficult because of atomic 
physics uncertainties. In addition, estimates of low optical depth in line ratios sensitive to resonance scattering 
effects also support a high-density emission scenario in the X-ray emitting regions of cTTS. 
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1. Introduction 

T Tauri stars are young pre-main sequence (PMS) late- 
type stars. "Classical" T Tauri stars (cTTS) are thought 
to still be surrounded by accretion disks as evidenced by 
IR and UV excess, while no signs for the presence of a disk 
are found in the so -called "weak line" T Tauri st ars (for 
details we refer to iFeieelson fe Montmerle|ll999|h X-ray 
emission from PMS stars is expected to be high because of 
their fast rotation if the emission is interpreted as scaled- 
up solar-type activity. However, for cTTS an additional 
source of X-ray emission through accretion is available; 
this additional X-ray production mechanism is expected 
to lead to significant differences in X-ray emission levels 
and variability, but in particular to differences in the spec- 
tral properties of the X-ray emission. 
High-resolution X-ray observations with the transmission 
and reflection gratings aboard Chandra and XMM-Newton 
have now been obtained for about two dozens of stars, 
but only for very few cTTS. The best data are usually 
obtained for the Ovu triplet located at 21.6 A ("r-line"), 
21.8 A ("i-line"), an d 22.1 A ("f-line") The f/i-line ratio 
is density-sensitive l|Gabriel fe Jordan] Il969|) . but in all 
cases no ratios below un ity are encountered for coronal 
sources l|Ness et al.ll2004) . In contrast, the available high- 
resolution spectra of cTTS show unusually low He-like f/i 
ratios in TWHya l|Kastner et al.ll2002HStelzer fc Schniittl 
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12004 and BP Tau l)Schmitt et al.ll2005|h The first obvi- 
ous conclusion was that the plasma in cTTS is at ex- 
tremely high densities suggesting its origin in an accretion 
shock rather than a "normal" magnetically active corona 
llKastner et alJl200^1 . However, f/i ratios could only be 
measured for Ovu and Neix, and the f/i ratios of those 
ions also depend on UV radiation fields if they are strong 
enough and located close to the origin of the X-ray emis- 
sion. Since the presumed accretion shock region is also 
expected to produce intense UV emission, the observed 
low f/i-ratios would then not contradict the accretion hy- 
pothesis, but need not necessarily imply high densities. 
The Chandra HE TGS spectrum of TW Hya has already 
been analyzed by iKastner et alJ l)2002j) with a variable 
abundance differential emission measure analysis and the 
identification and discussion of the anomalously low f/i- 
ratios in O vil and Ne IX. We address an alternative ap- 
proach to density d etermination using F exvii lines at 
17.05 A and 17. 10 A. iMauche et alJ l)200l|) were the first 
to use this ratio as density tracer in their study of 
the Chandra HETGS spectrum of the intermediate polar 
EX Hya, demonstrating that this ratio is considerably less 
sensitive to photoexcitation than He-like ions. Further, we 
investigate the effects of resonant line scattering which 
also depends sensitively on plasma density. 
The atomic physics of Fexvn, especially the 3c? — ► 2p 
(15 A range) and 3s — > 2p (17 A range) transitions, is quite 
complicated and extensive efforts have been spent with the 
conclusion that a number of indirect processes have to be 
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Fig. 1. HETGS spectrum and a parameterized best-fit model of TWHya (summed plus- and minus orders); line letter 
designations refer to the line parameters given in Tab. ^ 



considered apart from the standard collisional excitation 
(CE) theory. The inclusion of resonance excitation and 
inner-shell excitation from Fe xvi as well as radiative and 
dielectronic recombination from Fe xvin improved the sit- 
uation enormousl y, but re sidual discrepancies remain (for 
more details see |GvJ|2003) . In view of these difficulties we 
focus our analyses on the comparison of TWHya with a 
sample of stellar coronae, but refrain from any quanti- 
tative determination of densities. Comparison with theo- 
retical calculations is based on atomic data using APEC 
llSmith et al.1120011) 1 and the most recent calculations by 
ICJ (1200.^1 . 

2. Data Reduction and Analysis 

Chandra observations of TWHya (K7 Ve, d = 57 pc) were 
carried out with the High Energy Transmission Grating 
Spectrometer (HET GS; ObsID 5 48ksec, June 2000); 
details are given bv iKastner et all (|2002l) . We extracted 
count spectra using the Chandra Interactive Analysis of 
Observ ations (CIAO) software, but used our own tool 
CORA ijNess fc Wichmannl2 002') to determine line counts 
by fitting line templates and converted them to line fluxes 
using effective areas extracted with CIAO (using the tool 
fullgarf ). After examining the spectra of the plus and 
minus sides we use the sum (representing the effective 
area- weighted average) for our analysis, which can be done 
if no anomalies occur on either side (which is not the case 
to our knowledge). Any correction of line fluxes for absorp- 
tion is based on a value of Nn = 4 x 10 20 cm~ 2 , derived 
from a broad-band spectrum of TW Hya by Robrade et al. 
(2005). For comparison we extrated the HETGS spectra 
of various stellar coronae in exactly the same way to assess 
spectral differences between TWHya and purely coronal 
sources. 

3. Results 

For our analyses we measured fluxes for two sets of lines, 
five lines of Fexvn and the Lya and Ly/3 lines of H- 
like oxygen and neon. Fig. shows the HETGS spec- 
trum between 15-17.2 A illustrating the reliability of our 
line detections and flux measurements. The best fit above 
a source continuum of 20cts/A is also shown (FWHM 
0.016 A for all lines), and the derived counts and fluxes 
are listed in Table^together with the effective area values 

1 Version 1.3.1; available at http://cxc.harvard.edu/atomdb 



Table 1. Line flux measurements for TWHya. 



A [A] 


counts 


A [c ff m 

err 


ion 


transm 


M flux'"! 


3C 15.01 


35.5 ± 6.2 


16.0 


Fe xvn 


0.859 


7.15 ± 1.26 


3D 15.26 


32.3 ± 6.0 


23.2 


Fe xvn 


0.853 


4.45 ± 0.83 


[c] 16.00 


blend in O VIII 


Fe xvni 


0.835 


[1.15 ± 0.41] [c] 


[c] 16.07 


9.39 ± 3.4 


19.0 


Fe xvi 1 1 


0.834 


1.53 ± 0.55 


3F 16.78 


16.6 ± 4.4 


16.1 


Fe xvn 


0.816 


3.13 ± 0.83 


3G 17.05 


37.6 ± 6.4 


14.8 


Fe xvn 


0.809 


7.64 ± 1.30 


M2 17.10 


21.2 ± 4.8 


14.6 


Fe xvi i 


0.807 


4.36 ± 1.00 


10.23 


45.0 ± 7.3 


83.7 


Nex 




2.18 ± 0.35 


12.13 


200 ± 14.5 


43.1 


Nex 




15.8 ± 1.15 


16.00 


47.3 ± 7.1 


19.2 


O vin 




6.39 ± 0.97 


18.97 


121 ± 11.2 


8.87 


O vin 




29.8 ± 2.78 



la| Transmission efficiencies with A r H = 4 x 10 cm 

^lO -14 erg cm -2 s _1 ; corrected for JVh 

[c] ~ 75% of 16.07-A line blends with Ovm at 16.00 A 



used for conversion of line counts to fluxes and the trans- 
mission efficiencies for TVh = 4x 10 2 cm" 2 calculated from 
iBalucinska-Church fe McCammonl J 19921); we assume 
standa rd cosmic abundances from lAnders fe Grevessel 
(1989J). The Fexvn line fluxes are corrected for absorp- 
tion, but we did not correct the Lya and Ly/3 lines which 
will be investigated in detail below. In Table^we also list 
an Fe XVIII line at 16.07 A (2p 4 3s 2 P 5 / 2 to ground state) 
which we use to correct the Ovm Ly/3 line at 16.00 A 
to account for contamination by Fexvin at 16.004 A 
(2p 4 3s 2 P 3 / 2 to ground). The ratio of these two Fexvin 
lines varies slowly with temperature from 0.73 to 1.06 for 
logT = 6.3 - 7.3 (as predicted by APEC). Given the low 
temperature of TW Hya, the measured flux in O vin Ly/3 
was reduced by 75% of the flux in the 16.07-A line. We 
corrected the stellar Ovm Ly/3 line fluxes by interpolat- 
ing the slow temperature dependence using temperature 
estimates from the ratio of lines of Ovm at 18.97 A and 
Ovn at 21.6 A. The line fluxes in Tableware used to cal- 
culate the line flux ratios given in Tabled For the Fexvn 
line ratios the TVu-corrected fluxes were used while the 
Ly/3/Lya ratios are (not yet) corrected; the Ovm ratio is 
corrected for contamination with Fexvn, but blending in 
the Nex Lya line is ignored. 

We compare the line ratios of TW Hya with a large sam- 
ple of analogous measurements of stellar coronae in a va- 
riety of classes. We first focus on the ratio AA17. 10/17. 05 
as a function of density in Fig. TWHya is indicated 
by the light shaded area denoting the la uncertainty 
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Table 2. Line flux ratios in TWHya (notation as in 
Fig. El. 



AA 


flux ratio 


avH 


AA 


flux ratio 


avW 


M2/3G 


0.57 ± 0.16 


0.9 


3D/3F 


1.42 ± 0.46 


0.6 


3G/3F 


2.43 ± 0.76 


1.3 


3G/3C 


1.06 ± 0.26 


0.8 


M2/3F 


1.39 ± 0.48 


1.2 


3G/3D 


1.71 ± 0.43 


2.1 


3D/3C 


0.62 ± 0.16 


0.4 


3F/3C 


0.43 ± 0.19 


0.6 


0(/3/a) 


0.18 ± 0.03 


0.14 


Ne(/3/a 


) 0.13 ± 0.02 


0.14 



lions (ronoe lo s (T/K)-6.2..7.0) 
APEC; iog(T/K)=6.2 
APEC; iog(T/K)=7.0 
Gu (pri,. 

low-density limit by Gu (2003)*" 

TW Hyo » 



'"'error- weighted average of stellar measurements 

range from measurement errors. Stellar coronal densities 
are calcula t ed fro m Neix f/i ratios given for 18 stars by 
iNess et all l|2004l) . cleared of all the Fexix blending. A 
clear discrepancy between all stellar measurements and 
TWHya can be recognized. Theoretical predictions from 
APEC and Gu (priv. comm.; not yet including indirect 
processes), all in the temperature range logT = 6.2 — 7.0, 
are shown for two temperatures bracketing those temper- 
atures where the bulk of the Fexvn line formation is ex- 
pected to occur. While with APEC a quantitative deter- 
mination of a high density ~ 4 x 10 13 cm~ 3 is possible, 
the measurement of TWHya does not deviate from the 
low-density limit as predicted by Gu. The stellar measure- 
ments appear more consistent with the APEC-prediction, 
however, since all coronal sources are thought to be far 
hotter than log(T/K) = 6.0, also APEC underpredicts 
the measu red 1 7.10/ 17. 05-ratios. In contrast, the calcu- 
lations bv iGul (|2p03), including indirect processes (but 
providing only the low-density limit at different tempera- 
tures) , agree very well with the stellar measurements and 
an expansion of these calculations including all indirect 
processes as a function of density is likely to provide bet- 
ter estimates for TW Hya. Obviously, the observed Fe XVII 
17. 10/17. 05-ratio in TWHya is smaller than that typi- 
cally measured in coronal sources and it is larger than 
in EX Hya (cf. iMauche et alJl20o"ll) . suggesting a plasma 
density in TWHya larger than typically encountered in 
coronal sources, but smaller than 2 x 10 14 cm -3 as in- 
ferred for EX Hya. We also compare more ratios in Tabled 
with our stellar sample. In particular we found TWHya 
to be different in the ratios AA17.05/16.78 (3G/3F) and 
AA15. 26/15.01 (3D/3C), while in the other ratios no differ- 
ences between TWHya and the stellar sample are found. 
This suggests that the 17.05-A and the 15.26-A lines are 
anomalously enhanced in TWHya while those at 15.01 A, 
16.78 A, and 17.10 A show no peculiarities. 
A comparative analysis of the line ratios AA15. 26/15. 01 
an d (3D/3C) AA15.26 /16.78 (3D/3F) has been carried out 
by INess et all l)2003l) to investigate opacity effects. For 
both ratios TWHya shows significantly larger ratios pos- 
sibly indicating resonant line scattering effects. However, 
laboratory experiments suggest blen ding of the 15.26-A 
line with an Fexvi satellite line fe.g.. lBrown et ai]l200l|) . 
possibly explaining a trend of increasi ng AA15. 26/15. 01 - 
ratios with decreasing temperature ijNess et all l2003j) . 
Since TW Hya does have relatively low X-ray t emperature 
I Kastner et alJ l2002t IStelzer fc Schmittl |2004) , the same 




12 13 
iog(W,/cm" 3 ) 



Fig. 2. Density dependence of the Fexvn line ratio 
AA17.10/17.05 predicted by APEC and M.F. Gu compared 
with 18 coronal measurements and TWHya; plasma den- 
sities are derived from Neix f/i ratios. 



reason can account for the anomalously high flux in the 
15.26-A line. In view of this ambiguity we also studied the 
H-like Ly,3/Lya: line ratios of Ovin and Nex. In Fig. EJ 
we show a sample of stellar ratios of O vin with the blend- 
ing correction applied as described above, compared with 
the measurement of TW Hya and atom ic data pred i ctions ; 
the same line ratio has been used by iTesta et alJ l)2004t) 
finding II Peg and IMPeg to be anomalous. The TWHya 
measurement is marked with the hashed box bracketing 
the formal la uncertainties in the line ratio and in tem- 
perature. Theoretical predictions from APEC and Chianti 
agree quite well with each other and uncertainties in 
do not lead to different conclusions. However, we note 
in this context that the above blending correction is not 
straightforward, since it is predicted rather differently by 
Chianti and APEC (which we used for the correction). In 
view of the stellar measurements and the fact that the 
error bars do not include uncertainties from the blend- 
ing correction, this deviation from theory appears rather 
marginal. We carried out the same procedure for the Nex 
lines and found no deviation from the stellar measure- 
ments (Table 01. 

4. Discussion and Conclusions 

We identified and measured a new line flux ratio sensi- 
tive to density, which strongly supports earlier conclu- 
sions of high plasma densities in TWHya. The ratio of 
AA17. 10/17. 05 is a sensitive tracer of high densities with 
little contamination from UV radiation or temperature. 
Besides the extremely low f/i ratios in O vn and Ne ix we 
also found this ratio to be anomalously low compared to 
all stellar coronae. Unfortunately, quantitative constraints 
on density are still ambiguous because theoretical calcu- 
lations do not yet cover the full range of interactions be- 
tween the ground state and excited states. At any rate, the 
density of TW Hya appears to be higher than that of typi- 
cal stellar coronae, but lower than that of the intermediate 
polar EX Hya. Unambiguous X-ray density measurements 
as obtained here are also important for other cTTS. 
The 15.01-A and 15.26-A Fexvn lines provide a sensi- 
tive test for the effects of resonant line scattering. A com- 
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*LETG* measurements of 
Procyon, e Eri, and aCenB 



lag(T/K) 

Fig. 3. Theoretical predictions of Ly/3/Lya ratio as a 
function of temperature for Ovill and measurements of 
TWHya compared with stellar measurements (tempera- 
tures from Ovill/Ovn line ratios); shaded areas are pre- 
dictions for values of Ah < 4 x 10 20 cm~ 2 . 



pariso n of this ratio fo r TWHya with a sample of cool 
stars l|Ness et al.ll2003|) shows a large value for TWHya 
albeit with substantial error; similarly large values are 
found for EVLac and Prox Cen. However, as discussed by 
iNess et all l|2003l) the relative strength of the 15. 26- A line 
can also be explained by blending with low-temperature 
lines, leading us to conclude that the line ratio of the 
15.01-A and 15.26-A lines provides no unambiguous ev- 
idence for resonance scattering in TWHya. Next, we ex- 
amined the ratios of H-like Ly/3/Lya lines, which are also 
sensitive to resonant scattering (cf. Fig. |2J. While the 
Ly/3/Lya line ratio for TWHya is clearly above the the- 
oretical expectation, it does not differ significantly from 
those encountered in stellar coronae. Also, there is a 
nagging uncertainty about the blending correction with 
Fe xvin, so that any resonant scattering effects appear 
marginal. Since the situation is similar for Nex we con- 
clude that there is no clear evidence for any X-ray opti- 
cal depth effects in TW Hya and that the optical depth r 
should be below unity. 
The apparent absence of resonant scattering and the 




R.= 1.5OR 
d = 57.0pc 
!og(T D /K) = 6.40 
/=0.41 

" ' A i n) = -3.497 
!og(VEM(OVIII)) = 52.680 
iog(jS=1.161xl(r u *X«/WM/T )= -18.55 



log(N J on" 3 ) 

Fig. 4. Area vs. density N e from the requirement optical 
depth t < 1; area values for surface filling factors are 
indicated. 



value of the observed emission measure in a line can be 



combined as follows: We use the formula for the optica l 
depth t at line center derived by iBhatia fe Sabal l|200lf) 

t = 1.61 x 10~ 14 f\^Jlj^N ion L = f3N ion L where / denotes 
the oscillator strength and A wavelength, M atomic mass, 
Td temperature, A^ on ion density, and L path length. 
Next, the product Ni on L can be expressed in terms of 
the volume emission measure VEM= N^LA, where A" e 
denotes the electron density and A the area of the as- 
sumed cylindrical emission region. N e and A^n are related 
through N ion = 0.85N e A[ on , where A ion denotes the abun- 
dance of the considered ion relative to hydrogen, assumed 
to be 3.18 x 1(T 4 . We thus obtain VEM= T 08 ^|f n/r In 
Fig. 0| we plot (assuming r = 1 and 0.5) A as a function 
of N e ; note that this curve moves up for r- values below 
unity. Since the area shown is bounded above by (half) the 
stellar surface, we also plot a shaded area corresponding 
to surface filling factors of accretion hot spots between 0.5 
and 0.05. If these filling factors are indeed of the order of 
a few percent as usually assumed for cTTS, it is clear that 
densities in excess of 10 12 cm~ 2 are required to account for 
both the observed emission measure in TW Hya as well as 
the absence of any clear optical depth effects. 
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